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DETERMINATION OF ATMOSPHERIC STRUCTURE FUNCTION 


BY USING A SINGLE COHERENT DETECTOR 
I. INTRODUCTION 

Fried 1 has shown that atmospheric turbulence can impose severe restrictions 
on the design of telescopes for optical heterodyne detection systems. The 
atmospheric turbulence causes amplitude and phase fluctuations in the received 
signal which limit the effective telescope aperture size. This limitation is 
directly connected with the wave structure function and thus a knowledge of this 
structure function becomes important. 

Measurements of structure function have been proposed that use dual optical 
systems and dual coherent detectors 2 . Many times, however, the equipment 
available is a laser heterodyne receiver which contains a single beam optical 
system and a single coherent detector. In this paper we will describe and analyze 
a technique for using the single beam laser receiver to measure the wave struc- 
ture function. 

The measurement procedure is as follows: first, mask the telescope aperture 
with an opaque screen having two small holes, second, cover one hole and 
measure the output of the laser receiver and, finally, measure the output of the 
laser receiver with both holes uncovered. By repeating these two measurements 
for different hole separations, one can simply determine the wave structure 
function. 
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II. LASER RECEIVER DESCRIPTION 

A block diagram of the laser heterodyne receiver to be considered is shown in 
Figure 1. The receiver consists of a telescope used to collect the incident signal. 
The beam from the telescope is then beat with a local oscillator beam in a photo- 
mixer. The output of the photomixer is then passed through an I.F. amplifier and 
A.M. detected. 

The collection aperture of the telescope is located in the plane AA' as shown in 
Figure 1. The aperture is masked so that the incident signal can only enter 
through two small holes. The signal from these two holes is then focussed on 
the photomixer surface located in plane BB’. Concurrently, a local oscillator 
beam is directed onto the surface of the photomixer by means of a beam splitter. 

The photomixer which is shown in Figure 1 , gives an output current which is 
proportioned to the power incident upon its surface. More quantitatively, we 
have 

i p = 7 j' |E p ! 2 dA, y - constant (1) 

where E p is the total electric field on the surface and P is the surface area of 
the photomixer. The output of the photomixer is followed by an I.F. amplifier. 
This amplifier passes all frequencies on or near the beat frequency of the local 
oscillator and the incoming signal. Other frequencies are rejected. The output 
of the I.F. amplifier, i IF , is then A.M. detected. If we have 

ii F (t) - A(t) cos [A oj t + 0(t)] (2) 

where Ao> is the beat frequency and, A (t) and 0 (t) are slowly varying with 
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respect to Aw, then 


i d (t) ~ crA(t) a — constant (3) 

where i d is the output of the A.M. detector. 

III. THEORETICAL ANALYSIS OF RECEIVER RESPONSE 

We will now proceed with the theoretical analysis of the receiver input. The 
analysis will be kept as general as possible in this and the next section; then, in 
Section V , certain simplifying assumptions will be made to make the experimental 
measurement tractable. This procedure is useful since, if certain of the simpli- 
fying assumptions are not met exactly, one can use the general formulation to 
see how errors are introduced into the measurement. 


We will now calculate the field due to the incident signal in the plane of the 
photomixer. In Figure 2, we show a simplified drawing of the receiver optics 
without including the local oscillator subsection. The plane z = 0 is the aperture 
plane and corresponds to the plane AA' in Figure 1. It is an opaque plane with 
two arbitrarily shaped holes in it. The plane z = z p is the plane of the photo- 
mixer. The focusing and relaying optics are located between the aperture and 
photomixer planes, however, they are not shown for simplification purposes. 


If the opaque screen is illuminated by quasi-monochromatic light, E inc ( p , t) , 
then the field in the detector plane, E s (p, t) , is given by 

E S (P. t) = f h(£, p') E inc (£\ t) dp' (4) 

J Aj+ A 2 


where 


P ~ x Ax + ya y , p' ~ x'a x + y ' a y 


(5) 
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and h (p, p‘) is the two dimensional transfer junction from the z = 0 plane to the 
z = z P plane 3 . We have adopted this general formulation since it is not neces- 
sary at this point to specify h (p, p') in a more detailed manner. We should note 
that h (p, p') also depends upon z P , however, this dependence has not been ex- 
plicitly shown. The integration in Eq. (4) is an area integral over the two holes. 
This has been shown symbolically by writing A x + A 2 . 

In the absence of turbulence, we will assume that a plane wave is normally in- 
cident upon the aperture of the telescope, i.e., 

E. nc (p, t) = Ae^ u,t + ^ (6) 

where A and 4> are constants. The effect of turbulence will be taken into account 
by assuming that the amplitude and phase of the incident wave are distorted. The 
incident field becomes 

E inc (p, t) = A(£, t) e i [ ^ t+ ^- t ) ] (7) 

where the amplitude variation can be expressed in terms of the log-amplitude 
function ( (p, t) i.e., 

A s (£> t ) = A :,e f <£•*>. (8) 

Here A s is the r.m.s. amplitude of A s (P- t). We have chosen to normalize A s (p, t) 
to A s rather then to A since A s is more accessible to measurement. The time 
dependence of 1 and is due to the relatively slow changes in the gross proper- 
ties of the turbulence. The fact that these changes are slow compared with the 
time constants of the optical waves is what allows us to use Eq. (4) which was 
derived under monochromatic assumptions. 
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We will assume that the two holes in the opaque screen are small enough so that 
£(p, t) and 4> (/°. t) are approximately constant across each hole. Now plugging 
Eq. (7) into Eq. (4) and using the above assumptions, we obtain 

E $ (p, t) = Mj (p) e^ 1 W + j k“ st + + M 2 (p)e £ 2(t) + iUt+^>2( t )] (9) 

where 

Mi (p) = A s f h(p,p')dp' (10) 

2 " Jaj 

2 

and 


*i(t) = HP1' t), 
2 2 

(ii) 

<t > i (t) = 4> (pi , t) . 

(12) 


2 2 

The field, E p , in the plane of the photomixer consists of the sum of the signal 
field E s plus the oscillator field E 0 , i.e.. 


E P = E 0 +E s . (13) 

It will be assumed that the local oscillator field in the plane of the photomixer is 
given by 

E 0 (£, t) = A 0 (p) e j [ “ ot + *°<*' t >]. (14) 

Here we have included a time dependence in 4> 0 (p, t) to allow for the phase jitter 
of the oscillator signal. If we now use Eq. (13) in Eq. (1) , we have 

i p (t) = 7 j (E 0 + E s ) (E 0 + E s ) dp ( 15) 

i p (t) = 7 f E 0 E 0 dp + Re 2y f E s E 0 dp + 7 f E s E s dp (16) 

J D — Jd Jd _ 
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where a bar over a symbol means the complex conjugate of that symbol. For 
low noise heterodyne receivers, |eJ » |E s | . therefore 

i p (t) = y f E 0 E 0 dp + 2y Re fE s E 0 dp. (17) 

J D J D — 

Now plugging Eq. (9) and (14) in Eq. (17) gives 

i p ( t ) = y f A 0 2 O) dp ■ 

j D (18) 


Re | b i ( t ) e^ 1 + j + 4>i <*)! + Bj (t) + ^^ & ' t + 


<t> 2 (t)] 


where 


and 


Bi(t) = 
2 


2 j\ Mi (p) A 0 (p) dp 

J D 2 — 


A <x) - O) “ , 

s U 


(19) 


Next we write Bi(t) in polar form: 


Bi(t) = |Bi(t)| e^J (t) (20) 

2 . 2 


and use this expression in Eq. (18) to obtain 


i p ( t ) - y j A 0 2 (p) dp + | Bj ( t ) | e^ l(t) cos [A cut + 0j(t) + 0j(t)] 

J 0 (21) 

+ |B 2 (t ) | e e 2 (t) cos [Acut + 4> 2 ( t) + \p 2 (t )] 


Since the function Bi (t), £ 1 (t), and <$> j (t) vary slowly compared with Acu, Eq. (21) 
222 

can be viewed as a sum of a dc term plus two amplitude and phase modulated 
waves of frequency Aco. Upon passing the signal through the I.F., we just obtain 
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* IF (*■) 


B x ( t ) | cos [Awt + ^(t) + 

( 22 ) 

+ | B 2 ( t ) | e £ 2 (t) cos [Awt + 0 2 (t) + i/f 2 (t)]. 

Next the I.F. signal goes into an A.M. detector which has been described in the 
previous section by Eqs. (2) and (3). 

In order to apply Eqs. (2) and (3) , we combine the cosines in Eq. (22). By using 
this result, which is derived in Appendix A, we find 

i i F ( t ) = A (t) cos [Acdt + 9 (t)] (23) 

where 

A (t) = |b 2 | e^^)cos(^ 2+ ^ 2 )(24) 

The expression for 9 (t) is given in Appendix A but is not presented here since 
0(t) information is lost in A.M. detection. 

Now passing i IF (t) through the A.M. detector we find 

i d (t) = A(t). (24a) 

This signal is now amplified and recorded on tape for later processing. 

It should be recalled that the measurement of the total structure function will 
require not only a knowledge of the receiver output from the two hole aperture, 
but also the output when one hole is covered. This result can readily be obtained 
from Eq. (24) by setting B t = 0 or B 2 = 0. 
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IV. STATISTICAL CHARACTERIZATION OF INCIDENT WAVE 

The transmitted beam becomes statistical in nature after passing through the 
random or turbulent medium. As a result, the wave's log amplitude £ (p) and 
phase <P(p) become statistical quantities that must be described by probability 
distributions or characterized by statistical moments. In the past experimental 
and theoretical results have shown that under certain conditions £ and 4> obey 
independent Gaussian statistics 4 . This information can be used to calculate the 
wave structure function by using Eq. (24) . 


If we assume that the statistics of 4>{p) and £ (p) are locally homogenous and 
isotropic, then one can define the log-amplitude and phase structure functions as 

Dj> ( r ) = <J£ (Pi) - £ (£ 2 )] 2 )> (25) 

and 

v r > = <>(£i) - ^(p 2 )^y ( 26 ) 

In the above <u> is the ensemble average of u and r = \p x - p 2 \. The wave 
structure function D (r) is then defined as 

D (r ) = Dj?(r) + D^(r) (27) 


Now squaring Eq. (24) and taking the ensemble average, we have 


< A 2 > 


Lb, I 2 <e 2 ^!> + |b„| 2 <e 2 ^2> 


(28) 


+ 2 |Bj | B 2 1 +£ 2 > cos (<£j - <f> 2 + - <p 2 )^> 


By using Fried's 1 results: 
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(29) 


exp [21 (£>i)] = 

2 

and 

exp <^{f (pj) + £ (p 2 ) ± [<£(£i) ~ 0(p 2 )]})> = e~ ' A D(r ) (30) 

in Eq. (28) , we have 

<A 2 > = |Bj 2 + | B 2 1 2 + 2 (bJ | B 2 1 cos (i//j - 0 2 ) e -,,5D ( r ) (31) 

Therefore, we see that the mean square ensemble average of the A.M. detector 
output is directly related to the wave structure function. 

V. EXPERIMENTAL DETERMINATION OF STATISTICAL QUANTITIES 

We will now use the results of the previous two sections to devise a method of 
measurement of the wave structure function by using the output signal from the 

A. M. detector. 

As a first step, the output signal A (t) is sampled and <A 2 > is determined. We 
would now like to use Eq. (31) to determine D(r) from <A 2 >, however, the 
parameters IbJ, |b 2 | and 4> l - \p 2 are unknown. In order to simplify the 
measurement, we will require that the laser receiver be designed so Bj = B 2 = 

B. Under these conditions, Eq. (31) reduces to 

<A 2 > = 2 |b| 2 [1 + e-* D ( r >] (32) 

The measurement of D (r) can now be carried out as follows: first, cover one 
hole in the opaque screen and measure the output of the A.M. detector. This 
output will be denoted by A j (t). We now sample A x (t) and compute <A 2 >. By 
setting B t = 0 or B 2 = Oin Eq. (31) we obtain 
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2 


(33) 


B 


<Aj 2 > 


Next we measure the output of the A.M. detector with both holes open. Since we 
already know |B | from Eq. (33) we can readily find D (r) from Eq. (32). Now 
move the holes in the aperture mask so a new value of r is obtained and repeat 
the above procedure.* 


There are several optical arrangements which will make B 2 = B 2 . One par- 
ticular arrangement which is fairly insensitive to design inaccuracies is as 
follows: 


1. Adjust the local oscillator beam so that its amplitude and phase are 
constant across the photomixer surface P. 

2. Require that the signal passing through the holes in the aperture mask 
be focused in the plane of the photomixer. 

3. Make the small holes in the mask circular and of the same radius. 

4. Require that the diffraction pattern (Airy disc) from each hole be com- 
pletely within the photomixer surface P. 

5. Place the two circular holes so that they are an equal distance from the 
center of the telescope collection aperture. 

That conditions (1) - (5) produce B t = B 2 is verified in Appendix B. 


*The value of | B | will not have to be remeasured if the holes in the opaque screen are kept at a 
constant distance from the center of the telescope aperture. 
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In addition to determining the total structure function, the above procedure can 
also be used to determine the instantaneous phase difference if amplitude fluc- 
tuations are negligible. Let us assume £ j = £ 2 = 0 and Bj = B 2 . Then Eq. (24) 
becomes 

A(t) = y"2 |b i yr + cos A <f) (t) (34) 

where 

A0(t) = <p x (t) - cf> 2 (t) (35) 

Again, by covering a hole, the value of |b| is obtained and then when used in 
Eq. (34), the two hole measurement will yield the phase difference. 
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APPENDIX A 


We will now prove that two sinusoids with different amplitudes and phases can be 
combined into one sinusoid, i.e., 


cos (cot + /Sj) + a 2 cos ( cot + /3 2 ) - A cos (cot + 9) m 


To show that the above is true and find A and 6 in terms of and /3 1? we 

2 2 

rewrite the sinusoids on the left in their exponential form. Upon doing this, we 


obtain 




" e j(*>t+/3 2 ) + 

- ^ CQS ^ f 


\ 7 e jwt + \ y e~ }0)t 


- A cos (cot + 9) 


y - a x e } & 1 + a 2 e + ^ 2 


(a 1 cos ySj + a 2 cos /? 2 ) + j (a x sin J3 1 + a 2 sin f3 2 ) 


\y\ - [(a x cos j3 x + a 2 cos jB 2 ) 2 + ( a i + a 2 sin Z^) 2 ^ (5A) 


[a 2 cos 2 /3j + 2a 1 a 2 cos /3j cos fi 2 + a 2 2 cos 2 /? 


a 2 sin 2 /Sj + 2 a t a 2 sin /3j sin /? 2 + a 2 2 sin 2 /? 2 ]^ 


[a 2 + a 2 + 2 a 2 cos (J3 1 - (3 2 )] / 
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Since y = 


< 7 


tan~ 1 


a x sin /3 X + a 2 sin /3 2 
a x cos /3 1 + a 2 cos yS 2 


(BA) 


|y | e j <y , we have 

— |y | < ■y eJ " t + ^ 1 7 | e - -* < 7 e“ 1 wt - A cos (wt + 6 ) (9A) 

I7 1 cos (o;t + < 7) = A cos (wt + 0) (10A) 

A = I7I , <9 = <7 + 277 n , n = 0 , ± 1, ± 2, (HA) 
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APPENDIX B 




We would now like to show that Bj = B 2 , when conditions (1) - (5) of Section V 
are satisfied. By rewriting this condition and by using the definition of B j and B 2 
given in Eq. (19), we have 


K(£)A 0 (P>e^> dp = [ M 2 (p) A 0 ( / Q)e^^' t > d £ : (IB) 

J D J D 

If we now assume the oscillator beam amplitude and phase are constant across 
the photomixer surface, Eq. (IB) reduces to 


(£) dp - jj\t 2 (£) d £ 


(2B) 


Next we assume the system is focused, i.e., z p = f where f is the focal length of 
the system. Then the field in the photomixer plane is proportional to the Fourier 
transform of the aperture field. Therefore h (£.£'), given by Eq. (3) , is 


h(yO, P _ ' ) 


- i 2f P ■ (,P - 2 p ) 
ere - , CT - 


- j k f 


j X f 


(3B) 


and then using Eq. (3B) in Eq. (10), we obtain 


where 


M 2 x (£) 


cr e 


k 

2 1 P 2 


-i o 2 




P - 



(4B) 


If we now assume the holes are circular holes of radius a with centers located 
at p- and p 2 , Eq. (4B) becomes 
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where 


M, 1 (p) 


cr e 


•j 2f ^ 


| P a ( P' - PI 


) e 


+ j Ti* ■ £ 


dp' 


(5B) 


P a (P> 


1 , |p| < a 

0 , |p | > a 


(6B) 


and the integral is over the whole p' plane.* If we make the change of variable 

p" = p' - pi , Eq. (5B) becomes 
— 2 


Mi (p) - A s cr e 


. k r pf _ 1 

J f L 2 J? ' d 1 J 


(p") 


+ j 7_£> • P ' 


dp' 


(7B) 


Performing the two dimensional integral in cylindrical coordinates, we obtain 5 


where 


Mj 1 (P) = 


f (P) 


k 

+ j 7 p 


Pi 
“ 2 


f (p) 


k a f cr /kap\ -j 4? p 2 

— JilTf ) e 


(8B) 


(9B) 


and Jj (z) is the Bessel function of first order. Now using Eq. (8B) in Eq. (2B), 
we obtain 

f (p) e +J7 “ ‘- 1 dp = [ e +J7 ~ '- 2 f (p) dp (10B) 

J D J D 

If we now assume that the diffraction patterns of the circular holes, f (p) lie 
completely within P, then Eq. (10B) is unchanged if we extend the integration 
limits to infinity. Doing this (10B) becomes 

* Whenever integration limits are not specified, the limits should be taken are infinite. 


18 



(1 IB) 


J f (P) e+J 


k 

j£ ■ Pi 


dp - 

Now introducing cylindrical coordinates 

x = p cos 9 


f-** 


k 

J P * Pi 


f (p) dp 


y - p sin 6 


x i 

2 


yi 

2 


pi cos 6i 
2 2 


P l sin (9i 
2 2 


(12B) 


into Eq. (11B), we have 

• 03 27T . k 


n 2 7T “ 

e+) 

) 


P/Jj COS ( 0 ~ 0 l) 


f (yo) p dp d<9 

. r rv 

Jo Jo 


00 / ' 2?T - 1 - j f P P2 cos - ®2) 


(13B) 

f (p) p d p d<9. 


By using the integral representation for the Bessel function 5 


r 2 tt 

a) = e 

J o 


g+j acos (9 -y) jp 


2 77 Jo ( 


the angular integration in (13B) can be performed. The result is 


(14B) 


J* J 0 (j PP^j f (/0)pdp - J o ( J P P 2 J f O),°dp. (15B) 
We see if p x = p 2 , the two integrals will be equal and as a result 


Bj B 2 . 
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